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We present a Monte Carlo Renormalisation Group (MCRG) study of the SU(2) gauge theory with 
two Dirac fermions in the adjoint representation. Using the two-lattice matching technique we 
measure the running of the coupling and the anomalous mass dimension. We find slow running 
of the coupling, compatible with an infrared fixed point. Assuming this running is negligible 
we find a vanishing anomalous dimension, 7= —0.03(13), however without this assumption 
our uncertainty in the running of the coupling leads to a much larger range of allowed values, 
~0.6 < 7 < 0.6. We discuss the systematic errors affecting the current analysis and possible 
improvements. 
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1. Introduction 

Technicolor theories with fermions in higher representations of the gauge group can poten- 
tially provide a dynamical electroweak symmetry breaking mechanism without conflicting with 
electroweak precision data. Minimal Walking Technicolor is an example of such a theory, a SU(2) 
gauge theory with two Dirac fermions in the adjoint representation [1,2]. It is expected from per- 
turbation theory to be in or near to the conformal window, and lattice simulations have indeed found 
that the gauge coupling runs very slowly [3, 4, 5, 6, 7]. In order to be phenomenologically viable 
the theory must have a large anomalous mass dimension (7~ 1) [8, 9, 10]. A conjectured all-order 
beta function [11] predicts 7=11 /24 ~ 0.458, which is consistent with perturbative results in the 
MS-scheme up to four loops [12], and recent non-perturbative lattice results [3, 6, 13]. 

In this work we measure the running of the coupling and the anomalous mass dimension using 
the MCRG two-lattice matching method, a technique which has recently been used to investigate 
theories with many flavours of fermions in the fundamental representation [14, 15, 16]. 

2. Two-Lattice Matching Method 

We consider the evolution of all possible couplings of the system under Wilson Renormalisa- 
tion Group (RG) block transformations of scale s, where with each blocking step ultraviolet (UV) 
fluctuations ai^e integrated out, and irrelevant couplings flow towards their fixed point (FP) values. 
After a few such steps the only remaining flow is in the relevant couplings along the renormalised 
trajectory (RT). If we can identify two sets of couplings whose actions are matched after n and 
(« — 1) RG steps then their lattice coiTclation lengths must differ by a factor s. Since the physical 
conelation length should not be changed by the RG transform, this means that their lattice spac- 
ings, and hence inverse cutoffs A^' ~ a must also differ by a factor s. To identify such a pair of 
couplings, we need to show that after n and {n — \) RG steps respectively their actions are identical, 
which can be done by performing RG blocking on the gauge configurations and showing that the 
expectation values of all observables on these blocked gauge configurations agree. 

We use three 5 = 2 RG blocking transforms, labelled ORIG, HYP and HYP2 [14]. Each RG 
ti'ansform contains a free parameter a, and changing a changes the location of the FP, and how 
quickly the RT is approached in a given number of steps. Ideally it should be chosen such that 

• All operators predict the same matching coupling between {n,n — \) pairs for a given block- 
ing step n. (Deviations are a measure of the systematic eiTor from not being at exactly the 
same point along the RT.) 

• Consecutive blocking steps predict the same matching coupling, i.e. the coupling for which 
{n,n—\) pairs agree should be the same for all n. (Deviations show that the RT is still being 
approached in the irrelevant directions^) 

' If there existed other relevant couphngs that had not been precisely tuned to their critical values, these deviations 
could also indicate an unwanted flow in these relevant directions 
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3. Pure Gauge Results 

As an initial test of the method, matching in j3 between 32(16) and 16(8) lattices was per- 
formed using ~ 2000 SU(2) pure gauge configurations for each j3 . We matched in the plaquette, 
the three six-link loops, and three 8-link loops. An example of the matching of the plaquette is 
shown in Fig. la. The red, green and blue horizontal lines show the average plaquette on the 32^ 
lattice after 2, 3 and 4 ORIG blocking steps respectively, at j3 = 3.0, a = 0.57. The interpolated 
red, green and blue points show the average plaquette on the 16^ lattice after 1, 2 and 3 blocking 
steps respectively, as a function of j8', also at a = 0.57. The value of j3' where the two red lines 
intersect gives the matching coupling for n = 2, similarly the green and blue lines give the matching 
coupling for n = 3 and 4. 

This matching is repeated for each observable, and the spread of predicted matchings for each 
n gives a systematic error on the central matching value. The whole procedure is then repeated for 
various values of a, as shown in Fig. lb, to find an optimal value of a where subsequent RG steps 
predict the same matching value. The intersection of the last two blocking steps gives a central 
value for si,{^ = 3.0; 5 = 2) = jS — j8', while the range of couplings for which any of the blocking 
steps intersect within errors gives the uncertainty on this central value. Fig. 2a shows the resulting 
step scaling of the bare coupling Sb, determined using ORIG, HYP and HYP2 blocking on both 
32(16) and 16(8) lattices, along with the 1-loop and 2-loop perturbative predictions. 



SU(2) Pure Gauge; p=3.0: a=0.57: 32(16) plaquette matching ^-0' vs a, 32(16) lattice matching, p=3.0 




(a) Plaquette Matching (b) a-Optimisation 

Figure 1 : An example of the matching of the plaquette in jS for the pure gauge case using ORIG 
blocking on 32(16) lattices, and the a-optimisation including all seven matching observables. 



4. Discrete Beta Function Results 

Having confirmed that the method works for the pure gauge case, we now turn to the full 
Minimal Walking Technicolor theory of two Dirac fermions in the adjoint representation of SU(2). 
There are two couplings of interest in this theory, the gauge coupling and the mass. At an IRFP 
the gauge coupling is expected to be irrelevant, leaving the mass as the only relevant operator. In 
order to measure we first need to tune to the mass to its critical value, so that we are following 
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■ ■ ■ ■ l-loop 
2-loop 
O OORIG 32(16} 
O OORIG 16(8) 
O O HYP 32(16) 
O O HYP 16(8) 
O OHYP2 32(16) 
O O HYP2 16(8) 





(a) Bare step scaling Sj, for thie pure gauge tiieory 



(b) Bare step scaling Sj, for the full MWT theory 



Figure 2: The bare step scaling Sb for the pure gauge theory (consistent with perturbation theory), 
and the full MWT theory (consistent with a FP), using ORIG, HYP and HYP2 blocking. 



the running of the gauge coupling and not the mass, then the matching can be done in exactly the 
same way as for the pure gauge case. A fixed point would be indicated by a change of sign in the 
quantity S}, as the bare coupling is varied from weak to strong coupling. 

We use the HiRep [17] implementation of the Wilson plaquette gauge action with adjoint 
Wilson fermions and an RHMC algorithm with two pseudofermions. We then generated ~ 3000 
configurations on 16^ and S'* lattices for a range of jS values, each run at the critical bare mass. This 
allows two matching steps, after 2(1) and 3(2) steps on the 16'*(8'*) lattices respectively. 

The resulting measurement of Sb{^) is shown in Fig. 2b. This includes both the massless and 
some small mass {am ~ 0.005) runs; within errors Sb shows no mass dependence for these small 
masses. The ORIG matching values of Sh are clearly positive throughout, the HYP values are lower, 
and the HYP2 values are consistent with zero within error bai^s. There is no clear cross-over from 
positive to negative values of Sb for any of the blocking transforms, so while the data are consistent 
with a fixed point, they are not sufficiently precise to distinguish slow running from a fixed point. 

5. Anomalous Dimension Results 

Since the mass is the only revelant coupUng at the IRFP, we should in principle be able to 
match in the mass at arbitrary gauge couplings, given sufficient RG steps for the gauge coupling to 
flow to its FP value, and a mass small enough that despite the resulting flow in the mass the system 
remains close to the IRFP. In practice we only have a small number of RG steps, and because 
the beta function is small the coupling only flows slowly towards its FP value, so we set j8' = j8 
assuming the flow in the coupling due to one RG step is negligible. We match observables as for 
the pure gauge case, but instead of matching in j3 we match pairs of bare masses {amQ^a'm'o). 

We generated ~ 3000 configurations on 16^ and 8^ lattices, for a range of bare masses at each 
j8. This allows two matching steps, after 2(1) and 3(2) steps on the 16'*(8'*) lattices respectively. 

Because the bare mass is additively renormalised we convert the bare masses to PCAC masses. 
We measure the PCAC mass, am, as a function of bare mass, am^, for each j3 on the 16'* lattices. 
We then use this to convert the bare masses on both 8"* and 16"* lattices to PCAC masses, as the 
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measured PC AC masses on the 8^ lattices suffer from finite volume effects. Our previous result [18] 
for the anomalous mass dimension used PCAC masses measured on the 8^ lattices and hence 
contained a lai^ge finite volume effect, which has been removed in the present work. The anomalous 
mass dimension appears in the RG equation for the mass, 

d{am) 

^^^^ = -ymfl'M = -(1 + 7)am. (5.1) 

At an IRFP the anomalous mass dimension is a constant, so this can be integrated to give 

^^ = 2^+1 (5.2) 
am 

for a pair of matching masses {am,a'm'), from which a value for / can be extracted. We used four 
values of /3, j3 = 2.15,2.25,2.35,2.50, and the matching PCAC mass pairs using the HYP blocking 
tranform are shown in Fig. 3a. We also repeated the matching using ORIG and HYP2 blocking, 
and the agreement between the three blocking methods is very good. 

Different jS values seem to predict consistent values for the anomalous mass dimension, as 
shown in Fig. 3a, which uses all the /3 values and masses in the range 0.02 < am < 0.16. A 
combined fit to these data gives 7= —0.03(13), however significant systematic eiTors exist that 
have not been accounted for in this result, in particular the assumption that we can set j8' = j8. 



All p HYP matching tor 0.02<am<0.1 6 



Watching PCAC mass pairs for p=2.25 
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(a) Mass matching assuming j8 = j3' 



(b) Mass matcliing at /3 = 2.25 for a range of /3'. 



Figure 3: Matching masses in the range 0.02 < am < 0.16, with the assumption that j3' = j3, gives 
a vanishing anomalous mass dimension 7= —0.03(13). However allowing j3' to vary increases the 
allowed range of values to —0.6 < 7 < 0.6. 



6. Matching Observables 

In principle, for a given {P,am) there should be a unique matching set of couplings {P',a'm'), 
for which all blocked observables agree after n and (n — l) blocking steps respectively. In practice 
however, all of our observables are small Wilson loops, and as such are strongly correlated and have 
a very similar dependence on j3' and a'm'. This means that we can in fact find a "matching" a'm' for 
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a range of values of j3'. As an example the matching mass pairs for j8 = 2.25 and various values of 
j3' are shown in Fig. 3b. While we find that 5/, = /3 — j8' is compatible with zero, corresponding to 
setting j8 = /3', the error bars are relatively lai^ge, —0.08 ^ j8 — j3' < 0.16. From Fig. 3b we see that 
for j8 = 2.25 this region is approximately bounded by j3' = 2. 15 and j3' = 2.35, and encloses a large 
range of values for the anomalous mass dimension, —0.6 < 7 < 0.6. This range is representative of 
the errors in the anomalous mass dimension due to the uncertainty in the correct value of j8', and is 
the dominant source of systematic uncertainty in our results. 

To resolve this issue, more observables which are "orthogonal" to the Wilson loops, for ex- 
ample meson correlators or other fermionic observables such as the PCAC mass, would need to be 
included in the matching to determine which of these blocked configurations are actually matched, 
giving a unique matching set of couplings. 

7. Finite Volume 

The calculation of /requires us to compute the flow in the mass parameter for a conformal field 
theory (CFT) in the presence of a small mass deformation. Of course to be a true CFT the theory 
must be considered in infinite volume but unfortunately when we do simulations we have lattices 
of finite size, and so both the mass pai^ameter and the lattice size then determine any coiTclation 
length [19]. To extract the correct physics it is important to ensure that the correlation lengths we 
are measuring and matching are not being strongly influenced by the finite box size, m » 1 /L, 
however if we use too large a mass we will move the system a long way from any IRFP and the 
simple MCRG techniques we ai^e using will not apply. For our L=16 lattices to achieve for example 
mL > 2 we would need am > 0.12, which may be simply too large a mass to keep the system close 
to the fixed point. 

8. Conclusion 

We find a vanishing anomalous dimension and slow running of the coupling consistent with 
an IRFP, however our results suffer from considerable systematic eiTors. The most significant of 
these is the uncertainty in the location of the fixed point in the coupling, which in turn produces a 
large uncertainty in the value of the anomalous mass dimension. This is essentially the same issue 
that the Schrodinger Functional studies of this theory have suffered from. 

Assuming negligible running in the coupling we find 7= —0.03(13). Using our errors in the 
measurement of the discrete beta function to suggest a window for the possible flow in j3 produces 
a large uncertainty in the associated anomalous dimension —0.6 < 7 < 0.6. A small value of 7 is 
also likely to result from finite volume effects; if the beta function is small but non-zero then many 
blocking steps and hence large starting volumes will be needed to move the system away from the 
usual perturbative UV fixed point. 

Adding more matching observables, in particular fermionic ones such as meson conelators, 
should give a broader set of observables which will enable us to determine unambiguously matched 
actions. Using 32^ and 16^ lattices matching could be done after 1/2 and 2/3 blocking steps - 
the smallest blocked lattice being of size 4^, so that for example the meson correlator at each 
timeslice could be used as an observable. Going to 32^ lattices will also allow us to go to smaller 
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masses and thus be closer to the IRFP, as well as allowing us to match our current set of blocked 
observables after three blocking steps instead of two, which is a much more stringent matching 
condition that should help to constrain jS'. We will also be able to check for finite volume effects 
using the analysis of Ref. [16], and it may be possible to get a determination of the anomalous 
mass dimension using the stability matrix MCRG method. Other potential improvements include 
using an improved action to reduce ^(a) effects, and adding an adjoint plaquette term to the action 
to move the system away from a non-physical ultraviolet fixed point due to lattice artefacts, and 
allow us to go to stronger coupling [16]. 

However, even taking into account our currently large systematic errors, we find a small 
anomalous mass dimension which is clearly less than 1, and hence that, at least in its simplest 
form [20], the SU(2) gauge theory with two Dirac fermions in the adjoint representation is not a 
viable Walking Technicolor candidate. 
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